Emergent fermions and anyons in the Kitaev model 
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We study the gapped phase of the Kitaev model on the honeycomb lattice using perturbative con- 
tinuous unitary transformations. The effective low-energy Hamiltonian is found to be an extended 
toric code with interacting anyons. High-energy excitations are emerging free fermions which are 
composed of hardcore bosons with an attached string of spin operators. The excitation spectrum 
is mapped onto that of a single particle hopping on a square lattice in a magnetic field. We also 
illustrate how to compute correlation functions in this framework. The present approach yields 
analytical perturbative results in the thermodynamical limit without using the Majorana or the 
Jordan- Wigner fermionization initially proposed to solve this problem. 
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The study of elementary excitations in strongly corre- 
lated systems is a fascinating field of current research. As 
exemplified in the fractional quantum Hall effect, such 
excitations can be very different from the elementary 
constituents present in the system. In the same spirit, 
emergent fermions and gauge fields in purely boson/spin 
systems have attracted much attention recently [1]. The 
emergence of anyonic excitations in two-dimensional sys- 
tems has also triggered a tremendous amount of inter- 
est, especially since its relevance for topological quantum 
computation (see Ref. 2 for a recent review) has been 
pointed out by Kitaev in a seminal paper [3] introducing 
the celebrated toric code. 

More recently, Kitaev introduced a more realistic 
model [4] for which experiments using ultracold atoms 
or polar molecules have been proposed [5, 6]. This model 
is a two-dimensional spin-1/2 system on the honeycomb 
or brick- wall lattice, as illustrated in Fig. 1. It consists 
solely of Ising-like interactions but in different quantiza- 
tion axis. More precisely, the Hamiltonian reads 



H = 



E 



E 

a=x,y,z a — links 



Jn 



(1) 



where erf are the usual Pauli matrices at site i. Without 
loss of generality [4], in the following, we assume J a > 
for all a and J z > J x , Jy- 

Kitaev solved the model exactly by introducing Majo- 
rana fermions to represent the spin operators in an ex- 
tended Hilbcrt space. In this way, the Hamiltonian is 
reduced to free fermions in a static Z2 gauge field on the 
honeycomb lattice. The physical states are selected by 
a projection step which simply amounts to a selection 
rule on the parity of the number of fermions [7] , in agree- 
ment with the general conclusions of Ref. 8. The system 
exhibits a gapless phase for J x + J y > J z , with non- 
Abelian anyonic excitations arising when a magnetic field 
is switched on [4] . A contrario, for J x + J y < J z , the sys- 
tem is gapped and the low-energy effective Hamiltonian, 



at lowest nontrivial order in perturbation (J Xl J y <C J z ) 
and periodic boundary conditions, turns out to be ex- 
actly the toric code. Consequently, at this order, free 
Abelian anyons are present in the gapped phase. 
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FIG. 1: (color online). Mapping of the brick-wall lattice (a) 
to an effective square lattice (b) with unit basis vectors m 
and ri2. The numbering of the sites of a plaquette p is shown 
in both cases. 

Note also that an alternative treatment has been pro- 
posed [9, 10, 11] based on the Jordan- Wigner transfor- 
mation which transforms the spin system into a system of 
fermions with p-wave BCS pairing and a site-dependent 
chemical potential. 

In this Letter, we focus on the gapped phase and we 
derive the low- and high-energy effective theory at high- 
order in perturbation. This allows us to show that i) 
the low-energy effective Hamiltonian is an extended toric 
code Hamiltonian with still static but interacting Abelian 
anyons ; ii) the high-energy excitations arc free fermions, 
composed of a hardcore boson with an attached string of 
spin operators, hopping on a square lattice embedded in 
a magnetic field. 

Our approach has some advantages which compen- 
sate for its perturbative nature and its restriction to the 
gapped phase. It provides a unified picture of the emer- 
gence of high- and low-energy excitations without intro- 
ducing fermionic fields by hand as is done with Majorana 
or Jordan- Wigner fermionization. In addition, we work, 



2 



from the beginning, in the thcrmodynamical limit, and 
derive analytical results for non-translational-invariant 
low-energy states. 

Let us consider the model in the limit J x , J y <C J z . In 
the limiting case J x = J y = the model is a collection of 
isolated z-dimcrs. Each dimer has four possible configu- 
rations : two low-energy states {| ft), | II}} with energy 
— J z and two high-energy states {| Tl)> I It)} with energy 
J z . It is thus natural to interpret the change from a fer- 
romagnetic to an antiferromagnetic dimer configuration 
as the creation of a particle, with an energy cost that 
we set equal to 1 by chosing J z = 1/2. By construction, 
these particles are hardcore bosons hopping on the sites 
of an effective square lattice (see Fig. 1), together with 
an effective spin- 1/2 indicating which kind of (anti-)ferro 
dimer configuration is realized. Among the four possible 
mappings we choose the following 

irr> = ifro>,iu) = uM)), iu) = i^i),iit) = i^i), (2) 

where the left (right) spin is the one of the black (white) 
site of the dimer, and double arrows represent the state 
of the effective spin. Let us denote by b\ {b i ) the creation 
(annihilation) operator of a hardcore boson at site i (bold 
letters denote effective sites), and r" the Pauli matrices 
of the effective spin at the same site. With these nota- 
tions, the number of bosons in the system is Q = J^i 
and the Hamiltonian (1) can be rewritten as [7] 



H = - — + Q + T + T +2 + T_ 2 , 
where N is the number of z-dimers, 



(3) 



T ° = - E ( J * ^ + J y ^ + kc -) > ( 4 ) 

i 

t +2 = Y, ( J - u ' +ni + J y u ' +n2 ) = T -2' ( 5 ) 

i 

with hopping and pair creation operators 



fi+n-L _ jut T x 

L i ~ u i+m u i 'i+ni> 



= — iftt . b.r« 

1+n 2 % 2+T12 



„»+«■] 



— u i+n 1 u i '*+nii 



i+n 2 = ift t 



i+n 2 i i+n 2 



r?,(6) 



The vectors n\ and n,2 are shown in Fig. lb. Let us 
underline that [Wp,t?] = [Wp,u|] = for all (p,i,j) 
where the conserved plaquette operators [4] read (with 
notations given in Fig. 1) 

W p = alalalololal = (-1)^+^ r* r*. (8) 

At this stage, note that both the mapping (2) and the 
form of the Hamiltonian (3) are simply an alternative 
description of the problem which is always valid, even in 
a nonperturbative regime. The main difficulty resides in 
the fact that, now, one has to deal with hardcore bosons 
coupled to effective spin degrees of freedom. Of course, 



one could use a fermionization trick and solve the model 
directly as done by Kitaev [4]. However, the procedure 
used in the following can be applied to nonexactly solv- 
able models and further allows, in the present case, for a 
clear identification of the excitations. 

In general, a Hamiltonian of the form (3) cannot be 
diagonalized exactly. Here, following Kitaev, we choose 
to treat it perturbatively in the limit J Xl J y J z . Nev- 
ertheless, the Green's functions method used in Ref. 4 
turns out to be rather hard to implement at high or- 
der and/or high energy. Instead, we use an alternative 
approach based on continuous unitary transformations 
(CUTs) [12] whose perturbativc version [13, 14] is es- 
pecially well-suited to the problem at hand. Technical 
details will be given in a forthcoming publication [7]. 

The main idea is to transform the Hamiltonian (3) 
which does not conserve the number of bosons into an 
effective Hamiltonian H e s which satifies [H c s,Q] = 0. 
As explained in Ref. 15, i? c fr is a sum of fc-quasi-particle 
(QP) operators with k £ N. The k = 0, 1 contributions 
can be written as 



H. 



Oqp _ 



eff 



= Eo - Y 



, Pn w Pl w P2 



{Pl,...,p„} 



In 



■ r 



w v 



(10) 



(9) 



{ii>— >i„} 

Here {p\,p2, ■ ■ ■ , p n } denotes a set of n plaquettes and the 
sum, . . . ,j n } represents a sequence of n connected 
sites. The perturbativc aspect comes from the fact that 
the coefficients appearing in H e g are obtained as series 
expansions in J x and J y . We have computed the 0-QP 

amplitudes Eq and C Pl Pn up to order 10 and the 1- 

QP amplitudes fi and Dj ^ up to order 4. These 
expressions, being quite lengthy, will be given in a longer 
paper [7], as well as a discussion of (k ^ 2)-QP sectors. 

Let us first discuss the low-energy physics, i.e. the 
spectrum of H^ p . The most striking result is that the 
latter is only expressed in terms of the conserved quan- 
tities Wploqp = t^t^tIt^. Thus, for any (vortex) con- 
figuration of the W p 's which can take two values ±1, one 
readily gets the ground state energy of the correspond- 
ing sector. We wish to emphasize that at order 4 (lowest 
nontrivial order), as already discussed by Kitaev [4], the 
only nonvanishing contribution (apart from Eq) involves 
only single plaquette-terms and one recovers the toric 
code Hamiltonian with its Abelian anyons [3]. At order 
6, one gets the following corrections to these terms 



N 



</„ + J,. 



jy 



H — — T 2 T 2 
— 2 x V 



J X Jy 



j x j y 



j6 — i — 76 

(11) 



(12) 



but, more interestingly, one also has some two-plaquette 
terms : 



C 
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J X Jy 



-J 2 J 4 



(13) 
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Actually, when increasing the perturbation theory or- 
der, one generates higher and higher multi-plaquctte 
terms. Thus, the low-energy effective theory of the Ki- 
taev model turns out to be an interacting anyon theory 
whose eigenstates are those of the toric code. The single- 
vortex energy and the two vortices interaction energies 
of three configurations are given in Fig. 2 at order 10, for 
J x = J y = J. It should finally be noticed that although 




AE lv = J 4 + 8 J 6 + 75 J 8 + 784J 10 



AE 2v = 2AE lv - |j 6 - ±|?J 8 - ±fi j 



AE 2v = 2AE lv - if J 8 - 572,/ 10 



AE 2v = 2AE lv - ^ J 8 - 68 J 10 



FIG. 2: (color online). One-anyon and some two-anyon con- 
figurations (grey plaquettes) on a vortex-free background. 
AEi v (AZ?2v) is the energy cost for adding one vortex (two 
vortices) to the vortex-free state. We have set J x — J y — J. 

the vortices interact, they remain static since W p s are 
conserved quantities. 

In each sector defined by a configuration of the W p 's, 
we shall now see that the excitation spectrum is of 
fermionic nature. As proposed by Levin and Wen [8], 
the statistics can be directly determined from the hop- 
ping operators algebra. Let us consider the exchange 
process depicted in Fig. 3. The corresponding operator 
sequence is : 

itWM = -i, (14) 

or, equivalently, ijt^ = — tj^i*.. This latter identity 
shows that the quasi-particles made of a hardcore boson 
and an effective spin- 1/2 obey fermionic statistics. 

We shall now see that solving the excitation spec- 
trum is completely equivalent to solving a problem of 
free fermions in a magnetic field on the square lattice. 
Therefore, let us focus on the first-order perturbation 
theory for which H g F = Q + To and consider a single 
quasi-particle. Then, one can easily see that the mo- 
menta (Hlfi P — l) for all m are strictly equal to those 
of a pure hopping fermion Hamiltonian in a magnetic 
field provided the magnetic flux per plaquette mimicks 
the configuration of the Wp's. This is easily seen by not- 
ing that the product of the hopping operator around a 




FIG. 3: (color online). Illustration of the exchange of two 
particles discussed in the text, for j = i + n%, k — i + n\ and 
I = i — ri2- 

plaquette p is given by t\t l k t^tl = W p , with {i,j,k, 1} 
being any sequence of connected sites chosen among L, 
U, R and D (see Fig. 1). 

The Wp's being conserved quantities, from a purely 
spectral point of view one can replace the Pauli matrices 
in the hopping operators by pure numbers ±1 constrain- 
ing the fluxes per plaquette (in unit of the flux quantum) 
to be 4> p = for W p = +1 or <j> p = l/2 for W p = -1. At 
higher order, this mapping remains true (with hoppings 
from one site to any other site) but one needs to go be- 
yond this simple argument to prove the correspondence 
between both spectra [7]. Let us simply mention that 
this is done by building a basis of the 1-QP subspacc, 
which turns out to be made of states with one string 
of spin flips attached to one hard-core boson (yielding 
a fermion), as can be inferred from the form (6) of the 
hopping operators, and illustrated in Fig. 4. Note that 
the string fluctuates, as in the construction of anyons in 
[3]. 




FIG. 4: (color online). Pictorial representation of an emer- 
gent fermion as a composite object made of a hard-core boson 
(dot) and an attached string of spin-flips (thick line). Grey 
plaquettes are anyons, as in Fig. 2. 

To illustrate the mapping, let us consider the 1-QP 
vortex free state (W p = +1, Vp). The spectrum of -H e g P is 
made of a single band whose dispersion relation is given, 
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FIG. 5: (color online). Perturbative expansion at order 4, of 
the correlation function Cf 2 denned in the text, in the 0-QP 
subspace. 

at second order, by : 

E{ Iee (k x , k y ) = 1 - 2[J X cos(k x ) + J y cos(fc. y )] (15) 
+2[J X sm(k x ) + J y sin(fcy)] 2 , 

where the wave vector (k x ,k y ) £ [— tt,tt] x [— 7r,vr]. The 
gap in this sector is then obtained by minimizing Ef lee 
and one gets Af roc = 1 — 2(J X + J y ). Note that, in this 
case, the perturbative result at order 1 coincides with 
the nonperturbative result obtained by Kitaev [4] and 
one recovers the transition point at J x + J y = 1/2 = J z . 

Let us complete our analysis with a comparison to the 
Majorana (or Jordan- Wigner) fermion formalism. In the 
latter, one can compute the full spectrum exactly, by di- 
agonalizing, for each configuration of the W p 's with 2N 
spins, a 2N x 2N matrix. Within our approach, we have 
a N x N matrix to consider. Although we have seen 
that analytical (perturbative) results can be obtained 
for H®g P , even for non-translational-invariant Wp's, for 
Hlg P the problem is thus essentially as complicated as 
in the Majorana (or Jordan- Wigner) fermion formalism. 
In this case, the main interest of our method is thus to 
give a physical picture of the high-energy excitations as 
emerging fermions. 

Furthermore, the continuous unitary transformations 
allow for the calculation of correlation functions. In the 
0-QP subspace, as we show below, some of these can 
again be expressed as a series in W p operators, as was 
the case for the eigenvalues of H 3 p . Here, we focus on 
the correlation function C? z = (erf .erf G ) where i, • and 
i, o are the black and white sites of a given z-dimer i (see 
Fig. 1). Such a correlation function can be computed 
exactly for a translationnaly invariant configuration of 
Wp's, either by using the technique developped in [16] or, 
in a much simpler way, thanks to the Hcllman-Feynman 
theorem [7]. However, none of these two methods can be 
used effciently for a non-translational-invariant configu- 
ration of the WpS. 

As a first step, one has to translate C\ z into the ef- 
fective spin and boson degrees of freedom, which yields 
Cf 2 = ((-1) — )■ Th cn, the observable appearing in the 
angular brackets has to be transformed with the same 



unitary transformation as the one applied to the Hamil- 
tonian [7, 15]. As an illustration, we give the result in 
Fig. 5 at order 4 which is the lowest nontrivial order 
where the Wp's appear in the expansion. Note that we 
have used the same notation W for the operator and its 
expectation value on a 0-QP state. This result shows 
how the correlation function is modified by the presence 
of surrounding anyons. Of course, at higher order, anyons 
further apart would also contribute [7] as was already the 
case for the spectrum of H^g P . 

To conclude, we wish to emphasize that the method 
we have used here, based on a perturbative approach of 
the CUTs, can be applied to nonexactly solvable models 
as well. In addition, as we have shown, it is especially ef- 
ficient to compute the spectrum but also the expectation 
value of observables, even at high order in perturbation. 
Thus, we strongly hope that it constitutes a powerful tool 
to investigate, for instance, the non-Abelian anyons re- 
cently proposed by Yao and Kivclson for a time-reversal 
symmetry breaking version of the Kitaev model [7, 17]. 

We thank B. Doucot for fruitful and stimulating dis- 
cussions. 
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